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ABSTRACT 

The  important  feature*  of  nonequilibrium  flow  typical 
!  ■*  ( 

of  an  arc  tunnel  flow  system  are  outlined.  An  arc  heated  flow 
•yetem  is  analyzed  in  terms  of  its  subsystems  and  correspond- 
ing  local  flow  DamkShler  numbers.  Operating  ranges 
characterized  by  equilibrium  and  frozen  flow  limits  are 
specified.  The  effects  of  nonequilibrium  flow  on  test  results 
for  a  range  of  possible  operating  conditions  is  discussed. 

k 
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NOMENCLATURE 


A  Nozzle  section  area 

C  Atom  mass  fraction 

0^  Specific  heat  at  constant  pressure 

D  Diffusion  coefficient 

Dm  Damkohler  numbers 

h  Enthalpy 

h^  Dissociation  enthalpy,  ZC.h^ 

hQ  Total  enthalpy,  h^  +  h^ 

hD  Heat  of  recombination 

Constant  portion  of  recombination  rate  coefficient 

kw  sPe  cific  catalytic  surface  recombination  constant 

l  Characteristic  length  r*/tan  8 

Le  Lewis  number,  <r/Sc 

Mq  Molecular  weight  of  atomic  species 

m  Mass  flow  rate 

P  Pressure 

P  Total  pressure 

o 

q  Surface  convective  heat  transfer 

nc 

Surface  diffusion  heat  transfer 

q^.  Total  surface  heat  transfer  in  equilibrium  flow 

q  Total  surface  heat  transfer 
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NOMENCLATURE  (Cont.  ) 


Body  nose  radius 
Gas  constant  for  air 

Universal  gas  constant,  82.  06  (cc-atm/mole- *K) 

Nosale  section  radius 

Diffusion  velocity 

Entropy 

Particle  path 

Schmidt  number.  v/D 

Absolute  temperature 

Stream  velocity 

Compressibility  factor 

Degree  of  dissociation 

External  flow  velocity  gradient,  du^/dx 

Noasle  half  angle 

Thermal  conductivity 

Viscosity 

Kinematic  viscosity,  ^  /p 
Density 

Prandtl  number, 

Characteristic  time 

Correction  factor  for  catalytic  effects 
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NOMENCLATURE  (Concl. ) 


Subscript* 
c  Convection 

D  Atom  diffusion 

d  Dissociation 

E  Equilibrium 

e  Edge  of  boundary  layer 

f  Frozen 

i  ith  component  of  mixture 

N  Nitrogen 

o  Reservoir  state 

s  Stagnation  state 

w  Wall 

Superscripts 

* 


Nozzle  throat 
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I.  INTRODUCTION 


In  an  arc  tunnel  facility,  an  arc  heater  provides  high  enthalpy  gas  for 
long  periods  for  model  testing.  The  process  of  arc-heating  and  the  subse¬ 
quent  expansion  through  a  nozzle,  however,  introduce  a  nonequilibrium 
chemical  state  in  the  gas  in  the  test  section.  The  accuracy  with  which  this 
gas  state  is  known  often  has  a  direct  bearing  on  the  usefulness  of  the  arc 
tunnel  for  producing  test  results  of  quantitative  significance. 

Typically,  an  arc  tunnel  facility  consists  of:  (1)  an  arc  unit,  where 
energy  from  the  arc  discharge  is  introduced  into  the  gas;  (2)  a  plenum 
chamber,  where  the  gas  mixture  is  allowed  time  to  achieve  homogeneity; 
and  (3)  a  nozzle  through  which  the  prepared  fluid  from  the  plenum  is  dis¬ 
charged  to  form  a  well-defined  test  stream.  The  flow  around  the  model 
placed  in  the  test  stream  can  be  generally  characterized  by  a  shock  layer 
(when  the  test  stream  is  supersonic)  and  a  boundary  layer  through  which  the 
fluid  is  brought  in  contact  with  the  model  surface,  where  a  number  of  test 
measurements  are  made. 

The  chemical  state  of  a  reacting  flow  can  be  characterized  by  a 
Damkohler  number  Dm  which  is  a  measure  of  the  ratio  of  the  characteristic 
particle  resident  time  to  the  chemical  reaction  time.  Such  Damkchler 
numbers  can  be  defined  for  each  flow  region,  and  the  ranges  of  each  that 
correspond  to  equilibrium  and  frozen  limits  can  be  delineated. 

This  report  outlines  the  salient  features  of  the  nonequilibrium  flow 
typical  of  an  are  tunnel  flow  system  and  indicates  the  effects  of  the  non¬ 
equilibrium  flow  on  the  test  results.  Local  flow  states  in  separate  sub¬ 
systems  are  first  discussed  in  terms  of  their  local  Damkbhler  numbers. 

The  ranges  of  each  Damkdhler  number  that  correspond  to  equilibrium  and 
frozen  limits  are  delineated.  The  actual  determination  of  the  flow  states  in 
a  typical  arc  tunnel  is  demonstrated,  and  the  problems  of  establishing 
proper  operating  conditions  for  testing  are  discussed.  Other  problems 
that  arise  from  performance  of  a  test  and  interpretation  of  its  results  are 
also  d is  js sed. 


IL  SOURCE  FLOW  SYSTEM 


This  section  examines  the  chemical  states  of  air  passing  from  the  arc 

heater  to  the  plenum  and  through  the  expansion  nozzle  to  the  test  section. 

$ 

The  nominal  operating  ranges  of  typical  arc  units  are 

reservoir  stagnation  temperature  3.  500  -  8000*K 
reservoir  stagnation  pressure  0.  05  -  10  atm 

reservoir  total  enthalpy  3,  000  -  12.  000  Btu/lb 

The  chemical  properties  of  equilibrium  air  can  be  accurately  determined 
by  classical  thermodynamics  and  are  available  from  various  sources  (Refs. 

1  -  4).  Figu  re  1  shows  the  domains  of  major  chemical  activities  in  equili¬ 
brium  air.  Curves  showing  the  ranges  of  stagnation  pressure  and  tempera¬ 
ture  of  the  supply  gas  are  superimposed  on  these  domains,  indicating  the 
chemical  state  of  the  supply  gas.  It  is  clear  that  dissociation  of  oxygen  and 
nitrogen  are  the  important  chemical  activities  to  consider.  Energy  distribu¬ 
tion  in  the  various  degrees  of  freedom  is  given  in  Figure  2.  In  the  temperature 
range  3000  -  8000*K,  the  dissociated  species  have  high  energy.  Therefore  dis¬ 
sociation  and  recombination  kinetics  and  its  relaxation  effects  play  an  impor¬ 
tant  part  in  flow  determination.  Ionization  and  other  excited  electronic  states 
contribute  negligible  internal  energy  in  the  temperature  and  pressure  ranges 
of  interest.  For  aerodynamic  and  heat  transfer  simulations  these  effects  can 
be  ignored. 

A.  PLENUM 

The  plenum  allows  the  test  gas  to  achieve  mixing  uniformity  and  assures 
thermochemical  equilibrium  before  the  expansion  process.  The  detailed 
process  is  complicated  and  not  readily  describable.  To  provide  for  some 
gross  estimate  of  performance  we  can  establish  an  approximate  criterion  of 

* 

The  typical  operating  range  listed  here  represents  those  of  various  test 
facilities  currently  in  operation.  It  does  not  represent  those  of  facilities 
under  development  which  promise  to  achic  ve  higher  operating  ranges. 
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Figure  1.  Domain*  of  Major  Chemical  Activities  of  Equilibrium  Air 
Constituent*  and  Stagnation  Temperature  and  Pressures  Encountered 
as  Functions  of  Altitude  and  Velocity 
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Figure  2.  Energy  Distribution  in  Dissociation 
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chemical  time  necessary  for  equilibrium  by  correlating  the  available  results 
of  relaxation  studies  of  pertinent  chemical  processes  involved  (Refs.  7,  8, 
and  9).  An  upper  envelope  of  relaxation  times  required  for  air  constituents 
to  achieve  chemical  equilibrium  with  respect  to  dissociation-recombination, 
vibration,  rotation,  and  translation  was  obtained  empirically  based  on  the 
results  of  Refs.  7,  8,  and  9  and  is  shown  in  Figure  3.  This  envelope  can  be 
compared  with  the  estimated  residence  time  in  the  plenum.  When  residence 
time  falls  above  this  envelope  equilibrium  state  is  assured.  When  residence 
time  falls  below  the  envelope  some  modes  of  nonequilibrium  may  exist. 

Since  the  state  of  the  art  of  chemical  kinetics  of  air  constituents  is  6uch 
that  many  uncertainties  still  exist,  the  criteria  set  forth  above  can  be 
regarded  as  only  a  rough  guide,  to  be  corroborated  whenever  possible  by 
actual  measurements  of  gas  states.  A  Damkohler  number  for  the  plenum 


can  be  constructed 


<Pt) 


enveloi 


A 


where  residence  time  is  expressed  as  the  reciprocal  of  particle  velocity,  u, 
integrated  over  the  prrticle  path,  6.  The  chemical  state  of  the  gas  is 
indicated  by  the  following  limits: 


i  : 


Dm  j  <  1  equilibrium 
Drrij  »  1  nonequilibrium 


I 


The  correlation  is  such  that  only  the  equilibrium  limit  is  comparatively 
well  defined. 


Figure  3.  Approximate 

Air  i 


B.  NOZZLE  FLOW 


It  it  now  assumed  that  the  gas  reaches  thermochemical  equilibrium  in 
the  plenum  before  expansion  through  the  nozzle  proceeds.  The  rapid  accelera¬ 
tion  of  flow  in  the  nozzle  causes  high  rates  of  change  in  the  thermodynamic 
state  of  the  gas.  The  ability  of  the  chemical  state  to  keep  pace,  therefore, 
depends  on  the  relaxation  time  of  the  pertinent  chemical  processes  under 
local  flow  conditions.  In  general,  the  nozzle  flow  equations  must  be  solved 
simultaneously  with  chemical  rate  equations  to  obtain  an  accurate  account  of 
the  flow  variables.  A  number  of  machine  computing  programs  are  now  avail¬ 
able.  The  solutions  are  exact  in  the  sense  that  the  accuracy  of  these  results 
is  limited  only  by  the  correctness  of  chemical  reaction  rates  assumed  and 
by  simplifying  assumptions  that  are  us*»d  in  describing  the  nozzle  flow. 

For  complete  simulation  of  actual  atmospheric  flight  conditions,  it  is 
obviously  desirable  to  maintain  equilibrium  expansion  through  the  nozzle 
into  the  test  section.  In  most  shock  tunnel  applications,  efforts  have  been 
directed  toward  operations  under  equilibrium  flow  conditions.  This 
generally  requires  a  high  reservoir  pressure  and  a  relatively  low  reservoir 
temperature.  When  such  conditions  are  not  readily  obtainable  due  to 
operational  limitations  of  test  facilities,  the  use  of  nonequilibrium  test 
streams  must  be  considered.  Many  meaningful  test  results  can  be  obtained 
under  such  conditions,  provided  that  the  nonequilibrium  flow  states  can  be 
correctly  identified. 

lr,  high  enthalpy  arc  tunnels,  the  reservoir  pressure  level  is  limited 
by  arc-head  heat  transfer  limitations  or  available  electrical  power,  so  that 
nonequilibrium  nozzle  flow  is  a  primary  mode  of  operation. 

In  the  frozen  limit  of  a  nonequilibrium  flow,  the  chemical  states  of  the 
gas  are  fixed  and  the  determination  of  flow  is  therefore  greatly  simplified, 
provided  that  the  onset  of  the  frozen  states  can  be  correctly  specified.  Two 
simplified  methods  by  which  practical  calculation  of  frozen  flow  in  a  nozzle 
can  be  made  without  elaborate  machine  computation  are  discussed  below. 
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1.  SUDDEN  FREEZING  SCHEME 

Consider  a  diatomic  gas  wherein  dissociation-recombination  kinetics 
is  the  governing  mechanism  which  affects  the  departure  from  equilibrium. 

An  analysis  was  first  made,  with  particular  reference  to  nozzle  flow,  by 
Bray  (Ref.  9),  using  the  ideal  dissociating  gas  model  of  Lighthill  (Refs. 

10  and  11).  By  combining  the  one -dimensional  nozzle  flow  equation  and  the 
chemical  kinetic  relationships,  we  can  identify  a  Damktthler  number  for  the 
flow,  which  is  expressed  in  terms  of  local  equilibrium  thermodynamic  and 
flow  variables.  When  the  DamkBhler  number  is  larger  than  a  certain  value 
the  flow  can  be  considered  in  equilibrium.  On  the  other  hand,  when  the 
Damk&hler  number  drops  below  a  certain  value  chemical  kinetic  acitivities 
effectively  cease  and  the  chemical  constituents  of  the  gas  remain  constant 
(i.  e.  ,  "frozen").  In  a  nozzle  flow,  the  change  in  Damkohler  number  is 
drastic  enough  that  its  numerical  value,  corresponding  to  where  the  sudden 
freezing  occurs,  can  be  closely  bracketed,  giving  a  sudden  freezing  criterion. 
A  useful  simplified  scheme  for  nozzle  flow  solutions  has  been  devised  based 
on  such  criteria.  A  discrete  numerical  value  of  the  sudden  freezing 
Damkohler  number  was  obtained  by  Bray  from  correlation  of  the  result  of 
exact  solutions  (Ref.  9).  This  simplified  method  is  especially  useful  when 
knowledge  of  the  ultimate  chemical  state  of  the  gas,  but  not  of  its  inter¬ 
mediate  stages,  is  of  importance,  such  as  in  most  test  facility  applications. 
The  calculation  is  carried  out  by  assuming  equilibrium  flow  conditions  until 
the  local  DamkBhler  number  reaches  the  sudden  freezing  value.  Beyond  this 
point  calculation  proceeds  using  frozen  flow  conditions. 

Exact  calculations  of  high  enthalpy  air  flow  through  a  nozzle  of  a 
hypersonic  shock  tunnel  have  been  obtained  by  the  CAL  group  (Refs.  1Z,  13, 
and  14).  In  the  temperature  range  where  dissociation-recombination 
kinetics  is  predominant  (T  <  8000*K),  air  behaves  very  much  like  a  pure 
diatomic  gas  and  can  be  treated  ac  such  (Ref.  15).  Sudden  freezing  criteria, 
which  Bray  developed  for  a  diatomic  gas,  have  also  been  established  for 
oxygen  atom  recombination  in  air  (Ref.  12);  these  correlated  successfully 
with  the  exact  calculations  of  CAL. 
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The  freezing  criterion  makes  it  practical  to  obtain  results  with  a  desk 
calculator.  For  example,  Pearce  et  al.  (Ref.  5)  construct  a  map  of  estimated 
lrosen  atom  mass  fraction  in  terms  of  reservoir  states  of  air  for  a  specific 
nozzle  (1.  0  <  Pq<  1000;  4000  <  Tq  <  16,  000*K;  r*  =  0.  05  cm;  6  =  7.  5  deg) 
using  an  oxygen  recombination  rate  law.  This  map  summarizes  the  general 
trend  of  nonequilibrium  nozzle  flow  and  makes  it  convenient  to  fix  the 
reservoir  temperature  and  pressure  under  which  nonequilibrium  effects  are 
confined  to  tolerable  limits.  While  the  overall  trend  is  of  general  interest, 
the  specific  results  are  not  directly  applicable  to  the  arc  tunnel  of  present 
consideration  because  of  differences  in  pressure  range  and  nozzle  geometry. 
The  sudden  freezing  criterion  is,  in  essence,  an  empirical  scheme  based  on 
correlation  of  known  exact  results;  extrapolation  to  other  ranges  of  gas 
states  and  nozzle  geometries  must  therefore  be  considered  tentative  until 
iurther  checks  can  be  made  by  exact  results  in  the  corresponding  ranges  of 
gas  states  and  nozzle  configurations.  The  Aerospace  nozzle  flow  computing 
program  is  currently  being  adapted  to  obtain  solutions  directly  appropriate 
to  the  conditions  of  the  arc  facilities  under  development  in  this  laboratory. 

2.  FREEZING  POINT  ENTHALPY  CORRELATION 

Another  simple  method  of  predicting  the  freezing  level  of  dissociating 
atom  mass  fraction  due  to  oxygen  atoms  in  air  is  ou  lined  below.  Bray 
(Ref.  16)  has  made  further  study  of  the  freezing  criterion  and  noted,  from 
his  ideal  gas  model  with  a  single  chemical  reaction,  that  the  gas  enthalpy  at 
sudden  freezing  is  nearly  a  function  of  only  the  reservoir  entropy.  From 
the  results  of  previous  calculations,  he  has  constructed  on  a  Mollier  chart 
a  single  curve  on  which  the  freezing  states  for  all  reservoir  conditions  apper  r 
to  coalesce.  The  same  approach  was  extended  by  Bray  to  the  aii  expansion 
results  of  Hall  (Ref.  12)  for  a  predominantly  oxygen  atom  region  (Tq  =  5000 
and  6000 °K  and  Pq  =  100  and  1000  atm),  and  a  mean  freezing  curve,  which 
ha6  the  general  features  of  the  ideal  gas  case,  was  obtained.  In  Figure  4 
two  freezing  curves,  drawn  on  a  standard  Mollier  chart,  correspond  to 
normalized  throat  dimensions  bracketing  the  nozzle  dimension  range  of  our 
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interest.  These  curves  are  based  on  the  results  of  Ref.  16  and  are  extended 
to  the  region  of  present  interest.  A  few  cases  of  nozzle  calculation  by  sudden 
freezing  criteria  using  arc  tunnel  operating  conditions  correlate  with  the 
extended  freezing  curve  to  about  10  percent.  Since  the  dependence  on  the  i 
parameter  is  linear  (Ref.  16),  results  for  values  of  i  between  1  and  10  cm 
can  be  obtained  by  interpolation.  To  obtain  the  gas  freezing  state,  one 
needs  only  to  locate  the  stagnation  state  on  the  chart  and  drop  down  a  vertical 
(isentrope)  until  it  meets  the  freezing  curve.  The  intersection  gives  the 
freezing  state.  The  frozen  flow  downstream  of  t.ie  freezing  line  is  not  given 
on  this  chart.  For  a  given  fixed  freezing  curve  we  can  alter  the  downstream 
portion  of  the  Mollier  chart  so  that  it  is  appropriate  for  a  frozen  isentropic 
flow  (Ref.  5);  this  forms  a  combined  Mollier  -  frozen  flow  chart  for  use  in  a 
specific  nozzle  (Ref.  16).  Because  the  results  of  Refs.  12  and  16  apply  only 
for  the  case  of  oxygen  dissociation-recombination  kinetics,  the  plenum 
(starting)  conditions  of  Figure  4  are  limited  to  maximum  Z  values  in  the 
range  of  1.2,  below  which  nitrogen  dissociation  is  negligible. 

Williams  (Ref.  18)  derives  a  sudden  freezing  entnalpy  Damkohler 
number  Dm2  for  nozzle  flow  for  use  in  a  region  where  oxygen  dissociation- 
recombination  reactions  play  the  major  role  in  air  kinetics.  This  is  con¬ 
sistent  with  th»  discussion  in  the  preceding  paragraph,  and  the  same 
limitations  apply.  From  Ref.  18 


Dm2  = 


V*  '*** 


Mfcii  #fce 

a 


v 


where  the  oxygen  recombination  rate  coefficient  is  equal  to  k„(3000/T) 

jc  f  2  2  ^ 

with  k^  =  1.  16(10)  cm  -  *K  /mole  -cec.  Figure  5  shows  that  the 
enthalpy  at  sudden  freezing  is  largely  dependent  upon  Dm2  and  is  less 
sensitive  co  stagnation  temperature.  This  Damkohler  number  is  a  function 
of  easily  determined  plenum  and  throat  conditions  and  does  not  involve  the 
detailed  local  expansion  calculations  referred  to  in  Section  Bl. 
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Lordi  and  Mates  (Ref.  19)  discuss  calculations  of  chemical 
nonequilibrium  expansions  of  six  neutral  species,  N^i  O^,  NO,  N,  O,  and 
Ar,  which  are  vibrationally  and  electronically  excited  and  undergo  coupled 
chemical  reactions.  They  point  out  that  the  sudden  freezing  approximation 
is  not  usually  accurate  for  flows  with  coupled  reactions,  but  that  it  is 
possible  to  correlate  the  exit  frozen  enthalpy  with  plenum  entropy  for 
several  cases,  including  those  of  Ref.  14,  as  shown  in  Figure  6.  This 
correlation  from  Ref.  19  does  not  give  a  sudden  freezing  point,  but  rather 
tells  the  value  of  enthalpy  frozen  in  the  expansion  at  the  exit  of  large  area- 
ratio  nozzles  at  a  given  plenum  entropy.  Figure  6  is  used  to  locate  plenum 
conditions  and  frozen  exit  enthalpy  only  and  cannot  be  used  to  trace  the 
expansion  process.  Theoretical  cases  correlated  within  10  percent  in  the 
construction  of  Figure  6  involved  pressures  from  ten  to  several  thousand 
atm  and  temperatures  from  8000  to  15,  000*K.  As  is  shown  in  Section  IV, 
the  value  of  dissociation  enthalpy  frozen  into  the  flow  is  of  major  importance 
in  predicting  heat  transfer  rates  to  bodies  placed  in  the  flow. 


J o/R 


Figure  6.  Frozen  Enthalpy  at  Nozzle  Exit  Shown  on  a 
Mollier  Chart  for  Air 
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in.  FLOW  NEAR  THE  BODY 


This  section  examines  the  chemical  states  of  the  flow  about  a  test 
model  placed  in  a  supersonic  test  stream  provided  by  the  nozzle.  It  is 
assumed  that  the  Reynolds  number  of  the  nozzle  flow  in  the  test  section  js 
high  enough  that  continuum  flow  with  distinct  shock  layer  and  boundary  layer 
exist  in  the  neighborhood  of  a  body  (Re  =  p^u^R/^  »  25).  Although  in  the 
lower  range  of  the  Reynolds  number  some  solutions  for  viscous  shock  layers 
are  also  presently  available,  they  are  much  less  definitive  in  their  inter¬ 
pretation  and  will  not  be  included  in  this  discussion.  For  quantitative 
evaluation  of  test  results,  it  is  highly  desirable  to  operate  under  conditions 
where  the  flow  is  well  within  one  limiting  flow  regime  or  the  other 
(equilibrium  or  frozen).  This  is  true  for  nozzle  flow,  as  well  as  for  shock 
layer,  boundary  layer  flows,  and  surface  reactions.  In  the  following,  the 
flow  regimes  (equilibrium-frozen)  of  shock  layer,  boundary  layer  and  sur¬ 
face  are  outlined  separately  and  their  effects  in  the  determination  of  test 
results  are  noted. 

At  present,  arc  tunnel  operations  are  chiefly  geared  to  delivering  known 
heat  flux  to  models.  In  the  following  discussion,  the  effect  on  stagnation 
point  heat  transfer  will  be  considered. 

A.  SHOCK  LAYER 

Detailed  numerical  results  of  a  shock  layer  flow  in  front  of  a  blunt  body 
have  been  given  in  Ref.  20.  The  chemical  state  of  the  air  flow  in  the  shock 
layer  can  be  correlated  from  these  results  by  a  DamkShler  number.  Several 
similar  versions  of  the  Damkbhler  number  for  air  have  been  reported  (Refs. 

21  and  22).  The  one  given  below,  which  is  considered  most  current,  is  from 
Ref.  22. 


Dm  ^ 


exp(15.  8  -  60,  000/Tgf) 


(3) 


< 


where 


R  =  body  nose  radius,  ft 

=  free  stream  velocity,  ft /sec 
P  =  model  stagnation  pressure,  atm 
T  ^  =  frozen  stagnation  temperature,  *K 
X  =  universal  gas  constant,  82.  06  cc-atm/mole  *K 

The  practical  range  of  limiting  flow  regimes  has  been  indicated  as 
_2 

Dm^<  10  frozen  shock  layer 

Dm^  >10  equilibrium  shock  layer 

Among  flow  quantities  required  for  the  determination  of  heat  transfer 
at  the  stagnation  point  (see  Section  B)  all,  except  those  depending  on  surface 
conditions,  can  be  determined  from  shock  layer  analyses.  For  either  frozen 
or  equilibrium  flows,  they  are  readily  available  (Refs.  15,  23,  and  24).  It 
has  been  shown  that  flow  properties  in  the  shock  layer  are  not  sensitive  to 
the  dissociation  nonequilibrium  state  of  the  flow.  The  dissociation  level 
of  the  gas,  however,  needs  to  be  accurately  determined,  since  heat  transfer 
due  to  surface  recombination  in  a  frozen  boundary  layer  is  directly  propor¬ 
tional  to  the  energy  carried  by  atoms. 

B.  BOUNDARY  LAYER 

Similar  correlations  of  chemical  states  of  air  in  terms  of  DamkShler 
numbers  also  exist  for  the  boundary  layer  flows  (Refs.  15,  22,  25,  26,  27  and 
28).  A  comprehensive  correlation  derived  from  the  results  of  Refs.  15,  25, 
and  26  is  given  in  Ref.  27.  Such  results  are  for  noncatalytic  surface  condi¬ 
tions.  For  the  purpose  of  defining  flow  regimes  (equilibrium  or  frozen),  the 
correlation  is  valid  without  regard  to  the  surface  catalytic  activities.  The 
actual  heat  transfer  values,  of  course,  are  to  be  determined  in  conjunction 


where 


■"♦[£.(*•&)] 


(sec-  )  (using  hypersonic  approx.  ) 


T#e  =  stagnation  temperature  outside  boundary  layer,  *K 

h^  =  average  heat  of  recombination,  ZC^hR 

T  1 

_  2 

0^  dT  +  u^/2 


'‘■js'3 


with  the  surface  Damkdhler  number  to  be  discussed  in  Section  IV.  The 
boundary  layer  flow  Damkohler  number,  from  Ref.  27  is: 


=  mass  fraction  of  air  atoms  outside  boundary  layer 

=  mass  fraction  of  nitrogen  atoms  outside  boundary  layer 
1  e 

=  wall  temperature,  *K 

The  practical  range  of  limiting  flow  regions  is 
»-3 


Dm^  <10 
Dm^  >10* 


frozen  boundary  layer 
equilibrium  boundary  layer 
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IV.  SURFACE  CATALYTIC  REACTION  AND  HEAT  TRANSFER 


In  this  section  the  heat  transfer  to  the  model  of  catalytic  surfaces  in  a 
dissociated  stream  is  considered.  In  the  temperature  range  of  arc  jet  flow, 
the  energy  carried  by  the  atoms  in  dissociated  states  is  sizable.  The  amount 
of  dissociation  energy  released  to  the  surface,  in  turn,  depends  on  the  atom 
diffusion  flux  and  the  catalytic  efficiency  of  the  surface  to  promote  recombina¬ 
tion.  The  expression  for  heat  transfer  to  a  body  in  a  dissociated  stream  is 
first  given  as 


qw=qc  +  qD  (5) 

The  total  heat  flux  to  the  wall,  9W,  is  made  up  of  contribution  by  molecular 
convection,  qc>  and  that  by  atom  diffusion,  q^.  Expressions  iar  qc  and  q^ 
from  known  results  for  cooled  wall  blunt  body  stagnation  point  with  frozen 
boundary  layer  are  (Refs.  15,  23,  and  24). 


where 


I 


qc  =  (0.«)(2(V„Pt,)‘/27^/5  hefffw«h0>  (6) 

qD  =  <0.47)(2ftie,Pe,)1/2Sc"*/3V 


[1  +TT1'1-  »<*<  1  (0<^S<») 

(o.47)  <zi*e>pei> 1/2  se2/3pw_1 

specific  surface  catalytic  recombination  coefficient 

•*p 

rcihR.’  * =  ~r  •  Sc  =  pb 
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The  factor  0  it  derived  from  the  simultaneous  solutions  of  momentum  and 
species  -equations  in  a  irosen  boundary  layer  flow  of  dissociated  gases  with 
catalytic  surface  recombination  {Ref.  24).  It  can  be  considered  a  surface 
catalycity  parameter  whose  limiting  values  correspond  to  that  of  the  surface 
catalytic  activities  in  terms  of  k^,  as  noted  in  the  above  expression.  Con¬ 
sistent  with  the  above  discussion,  the  net  effect  of  surface  catalytic  reaction 
on  heat  transfer  can  also  be  discussed  in  terms  of  a  surface  Damkohler 
number,  which  is  related  to  0  as  follows: 

k 

Dm5  =  -gr-  (0<Dm5<  oo)  (7) 

As  noted  in  the  general  discussion  of  Dam k Shier  number  outlined  in  the 
Introduction,  surface  catalycity  alone  is  not  sufficient  to  establish  a 
criterion.  The  overall  effect  depends  as  well  on  the  flow  characteristics  in 
terms  of  quantities  in  the  expression  of  5,  the  characteristic  diffusion 
velocity  across  the  boundary  layer.  The  range  of  Dm,,  attainable  in  practical 
conditions  has  been  discussed  in  Refs.  22  and  24.  The  specific  surface 
catalytic  reaction  coefficient  is  a  material  property  which,  for  each  sur¬ 
face  material,  depends  on  the  gaseous  reactant  and  surface  temperature. 
Some  typical  values,  obtained  experimentally,  are  available  in  Ref.  24. 

A  few  approximate  analyses  of  nonequilibrium  boundary  layers  with 
arbitrary  surface  catalytic  activities  are  available  (Refs.  29  and  30).  They 
are  useful  in  assessing  the  relative  effects  of  the  boundary  layer  nonequili¬ 
brium  and  the  surface  catalytic  activity  when  the  limiting  conditions 
(equilibrium  or  frosen  boundary  layer,  fully  catalytic  or  non-catalytic  sur¬ 
face)  are  not  achieved  because  of  the  limitations  imposed  either  by  the  test 
facility  or  the  choice  of  surface  material.  Results  of  Ref.  29  are  given  here 
as  Figure  7.  Values  of  q^  -  qc/qE  -  qc  are  presented  in  terms  of  Dm^  and 
Dm,.  where  q^  is  the  actual  total  heat  transfer  to  the  cooled  blunt  body 
stagnation  point  and 

qE  =  0.47(2ftip)1/2a*2/3(ho) 


(8) 


•  ^ 

i*  the  heat  transfer  for  equilibrium  flow.  The  boundary  layer  correlation  for 


frozen  flow  (D^  <  10~*)  is  included  in  Figure  7,  and  by  use  -of  Eqs.  (5)  and 
(6)  can  be  shown  to  be 


V.  DETERMINATION  OF  FLOW  STATES  IN  A  TYPICAL  ARC  TUNNEL  AND 

EXPERIMENTAL  CONSIDERATIONS 


Arc  tunnel  flows  are  often  required  to  operate  under  nonequilibrium 
states  in  some  part  of  the  flow  system.  Correct  determination  of  the  flow 
in  each  of  the  subsystems,  as  outlined  above,  is  therefore  essential  for 
meaningful  interpretation  of  the  test  results  obtained.  The  typical  operating 
mode  of  an  arc  tunnel  entails  an  equilibrium  plenum  and  a  frosen  nozzle 
flow.  Options  are  open  in  principle  for  the  flow  about  the  model.  The  shock 
layer  and  boundary  layer  each  can  be  separately  in  equilibrium  or  frosen 
states.  But  in  practice,  the  prevailing  flow  conditions  are  such  that,  in 
most  instances,  both  shock  layer  and  boundary  layer  are  likely  to  exist  in 
equilibrium  or  frozen  states  concurrently.  In  the  equilibrium  case,  the 
fluid  state  is  determinate  without  regard  to  the  flow  history  prior  to  its 
passage  to  the  shock  layer.  If  at  the  same  time  the  surface  temperature 
is  sufficiently  low,  surface  catalycity  is  also  irrelevant  to  the  result  of 
heat  transfer,  since  equilibrium  atom  concentration  is  insignificantly  low  at 
low  temperatures,  so  that  few  atoms  will  reach  the  surface  before  recombi¬ 
ning  in  the  boundary  layer.  In  the  frozen  case,  however,  the  energy  carried 
by  the  dissociated  species  is  delivered  directly  to  the  surface  and  released 
through  catalytic  recombination  at  the  surface.  Although  the  differences  in 
fluid  properties  as  a  result  of  differences  in  flow  states  are  considered  of 
minor  importance,  the  surface  catalytic  activity  plays  a  significant  role  in 
the  determination  of  surface  heat  transfer. 

In  Sections  A  through  D,  criteria  discussed  above  are  used  to  determine 
the  representative  operating  conditions  that  exist  in  the  200-kw  supersonic 
arc  facility  presently  in  use  in  this  laboratory.  Experimental  considerations 
regarding  heat  transfer  measurements  under  such  conditions  are  also 
discussed. 
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The  reservoir  density  is  estimated  from  equilibrium  conditions  to  be 


The  velocity  is  estimated  from  the  measured  mass  flow  rate  and  the  sonic 

throat  conditions.  A  conservative  estimate  of  residence  time  in  the  plenum 

is  made  by  taking  the  velocity  as  given  by  the  one  dimensional  flow  through 

the  plenum  in  the  axial  direction.  Without  making  allowance  for  mixing, 

which  undoubtedly  is  taking  place,  this  represents  the  shortest  mean  particle 

path  for  the  flow.  The  measured  mass  flow  rate  is  1.  86  g/sec,  from  which 

resident  time  is  estimated  as  about  300  psec,  which  gives 

(Pt)  .  .  a  0.  6  X  10"*.  This  point  lie^s  above  the  envelope  of  Figure  3, 
resident 

indicating  tl.at  equilibrium  is  achieved  in  the  plenum.  At  lower  reservoir 
pressure  and  temperature,  conditions  are  increasingly  marginal.  Since  it 
is  difficult  to  obtain  sensible  estimates  of  resident  times  under  actual  plenum 
flow  conditions,  development  of  some  experimental  sensor  may  be  desireable 
to  provide  for  a  more  direct  method  of  obtaining  the  information  required. 


B. 


NOZZLE  FLOW 


Calculations  based  on  the  sudden  freezing  scheme  outlined  above 
indicate  that,  under  the  operating  range  of  the  nozzle  in  this  laborator 
freezing  occurs  in  most  instances  very  near  the  sonic  nozzle  throat.  I'hsre- 
fore  a  substantial  amount  of  rgy  is  carried  by  the  frozen  species.  The 
magnitude  of  this  problem  is  illustrated  for  a  supersonic  nozzle  in  Figure  8. 


V 


*  4. 


MC-jf. 
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The  chemical  reactions  are  assumed  to  freeze  at  the  throat  and  Le  =  1.  It 
is  seen  that  for  models  having  the  maximum  spread  of  effective  catalytic 
efficiency  (*=  0  to  *  =  1)  large  differences  in  heat  transfer  are  possible.  A 
model  with  *=  1  would  be  expected  to  receive  over  3  times  the  heat  flux  to  a 
+  =  0  model  at  the  highest  enthalpy  illustrated. 

C.  SHOCK  LAYER  AND  BOUNDARY  LAYER 

The  chemical  states  of  shock  layer  and  boundary  layer  are  near  enough 

that  their  operating  conditions  are  governed  by  nearly  the  same  ranges  of 

flow  parameters.  Only  boundary  layer  calculations  are  given  here.  Figure  9 

shows  approximate  values  of  the  boundary  layer  flow  Damkbhler  number 

Dm^  for  a  M  =  5  flow  over  hemispherical  models  for  a  range  of  plenum 

enthalpy  and  pressure.  Even  for  the  maximum  plenum  pressure  of  0.  5  atm, 

as  found  with  the  200-kw,  M  =  5  arc  in  this  laboratory,  little  gas  phase 

recombination  can  be  expected  for  models  in  the  range  of  1-in.  radius,  as 

_2 

Dm^  much  above  10  is  not  possible  (see  Figure  7).  As  Dm^  ~R,  it  is 
clear  that  probes  of  a  smaller  radius  could  only  decrease  Dm^  *o  give  more 
assurance  of  freezing.  On  the  otl  nr  hand,  much  larger  probes  are  limited 
by  nozzle  exit  diameter,  so  that  the  conditions  of  frozen  shock  and  boundary 
layers  are  rigidly  fixed  for  this  particular  arc  tunnel. 

D.  SURFACE  CATALYCITY  AND  HEAT  TRANSFER 

The  plot  shown  in  Figure  10  is  convenient  for  estimating  possible  heat 
transfer  differences  over  the  range  of  and  (h^/hQ)  expected  in  any  facility. 
This  plot  illustrates  the  total  heat  transfer  vo  a  model  assuming  no  gas  phase 
recombination.  Figure  10  incorporates  the  simplifying  assumption  that  Le  =  1 
in  Eqs.  (5)  and  (6),  which  can  be  shown  to  give 


IwMe  =  1+  <*-  l)j- 


It  is  of  interest  to  compute  limiting  values  of  ?  for  the  Mach  5  arc  tunnel  in 
this  laboratory.  Values  of  k^  quoted  in  the  literature  range  from  1.  0  cm/sec 
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Figure  9.  Boundary  Layer  Damktthler  Number  for  a  Range  of  Plenum  Enthalpy 


Figure  10.  Heat  Flux  Ratio  as  a  Function  of  h^/1^  and  ? 


for  noncatalytic  surfaces  to  over  10  cm/sec  for  the  best  catalyst.  Very 

small -diameter  hemispherical  probes  (R»  0.  01  in. )  would  yield  values  of  S 

3 

in  the  range  of  6  X  10  cm/sec  maximum,  while  minimum  3  values  for 

2 

usable  large,  flat-faced  models  (R  *  1.  0  in. )  would  be  near  3X10.  Thus, 
the  combination  of  a  small  diameter  probe  with  a  noncatalytic  surface  gives 
Dg  <  10"^,  or,  effectively,  *  *  0.  A  large  catalytic  model  gives  Dg  >  10, 
or,  effectively,  *  ■  1.  This  illustrates  the  range  of  conditions  possible  with 


the  present  facility.  Many  references  quote  values  of  k^  near  10^  cm/sec 
for  many  eigineering  materials,  so  that  with  intermediate* sized  probes  it 
is  not  unreasonable  to  at  times  expect  values  of  ?  in  the  0.  50  range. 

The  above  example,  used  in  conjunction  with  Figure  10  for  values  of 
(h^/ho)  up  to  0.  70,  serves  to  point  up  the  wide  range  of  heat  transfer  varia¬ 
tions  to  be  expected.  As  mentioned  in  Section  HI,  it  is  often  advantageous  to 
operate  under  limiting  conditions  with  ?  -  0  or  1.  However,  one  is  seldom 
completely  assured  of  this  condition,  since  k^  values  reported  for  a  single 
material  vary  over  a  range  of  three  orders  of  magnitude,  depending  on  the 
facility  and  test  techniques  employed.  There  appears  to  be  a  need  for  develop¬ 
ment  of  a  "self  calibrating"  probe  technique  along  the  lines  of  that  discussed 
in  Ref.  31  that  can,  in  fact,  be  used  independently  of  kw  values  determined 
elsewhere.  Reference  32  presents  an  experimental  study  wherein  the  sur¬ 
face  properties  of  arc  tunnel  probes  were  calibrated  independently  in  a  glow 
discharge  shock  tube  facility  of  known  dissociation  level  before  and  after 
arc  testing.  This  technique  was  found  satisfactory  for  determining  arc  jet 
oxygen  dissociation  levels. 

The  importance  of  the  parameters  (h^)  and  V  becomes  increasingly 
clear  when  one  is  faced  with  the  necessity  of  calibrating  the  stream  for  model 
heat  flux  and  enthalpy  level  using  a  calorimeter  of  a  catalytic  efficiency 
different  from  that  of  the  actual  test  model.  This  situation  can  become  even 
more  involved  if  it  is  desired  to  perform  meaningful  testo  of  an  identical 
nature  in  several  arc  tunnels.  In  this  case  one  faces  the  strong  possibility 
that  each  tunnel  may  operate  with  a  different  atom  concentration  at  the  same 
heat  flux  to  a  common  calorimeter.  Different  heat  transfer  rates  then  appear 
during  final  tests  in  each  tunnel  because  models  with  different  effective 
catalytic  efficiencies  (?)  are  operated. 
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VL  CONCLUSIONS 


It  has  been  shown  that,  because  of  physical  limitations,  arc  tunnels  are 
often  required  to  operate  with  nonequilibrium  flows.  Unlike  equilibrium  flow, 
for  which  chemical  states  in  the  flow  are  uniquely  determined  by  local  thermo¬ 
dynamic  parameters,  nonequilihrium  flow  necessitates  monitoring  the 
chemical  states  throughout  the  entire  system.  The  monitoring  is  facilitated 
by  the  use  of  Damkohler  numbers,  in  terms  of  which  limiting  flow  regimes 
can  be  defined.  In  the  frozen  limits,  flow  quantities  are  again  well  defined. 
Interpretation  of  test  results  obtained  therein  is  therefore  simplified. 

In  surface  heat  transfer  measurements  made  in  a  frozen  nozzle  flow, 
both  the  frozen  dissociation  level  and  the  surface  effective  catalycity  are 
important  factors  whose  effects  on  heat  transfer  has  been  demonstrated  by 
sample  calculations  based  on  the  operating  condition  of  the  arc  facilities  of 
this  laboratory.  Reliable  values  of  these  parameters  are  therefore  essential 
for  meaningful  interpretation  of  the  measured  heat  transfer,  especially  when 
measurements  are  made  in  different  facilities. 
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